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ABSTRACT 



We correlate the positions of 13,240 Brightest Cluster Galaxies (BCGs) with 
0.1 < z < 0.3 from the maxBCG catalog with radio sources from the FIRST 
survey to study the sizes and distributions of radio AGN in galaxy clusters. We 
find that 19.7% of our BCGs are associated with FIRST sources, and this fraction 
depends on the stellar mass of the BCG, and to a lesser extent on the richness 
of the parent cluster (in the sense of increasing radio loudness with increasing 
mass). The intrinsic size of the radio emission associated with the BCGs peaks 
at 55kpc, with a tail extending to 200 kpc. The radio power of the extended 
sources places them on the divide between FR I and FR II type sources, while 
sources compact in the radio tend to be somewhat less radio- luminous. We also 
detect an excess of radio sources associated with the cluster, instead of with the 
BCG itself, extending out to ~ 1.4 Mpc. 

Subject headings: galaxies: active — galaxies: clusters: general 



1. Introduction 



Recent simulations suggest that the colors of massive galaxies in the local Universe can 



only be explained if AGN feedback quenches star formation in the host galaxy (ICroton et al. 



2006). The link between supermassive black hole and galaxy formation (IHaehnelt fc Kauffmann 
20001 ) may depend critically on AGN feedback, via coupling of the mec hanical power (via 
wind s or jets) of AGN to the baryonic compo nent of forming gal axies ( Rawlings Sz Jarvisl 
2004f) . and radio jets are known to both quench ( Springel et al. 2005 ) and enhance ( Croft et al. 



20061 ) star formation. Mergers, galaxy harassment, and other processes invoked in AGN trig- 
gering, are common in clusters. 
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Studies of the AGN and galaxy populations, kinematics, and substructures in clusters 
provide insight into how galaxies, AGN and clusters of galaxies form and evolve. 

Clusters also contain some of the most luminous galaxies in the local Universe, the 
Brightest Cluster Galaxies (BCGs), which are more likely to be radio-loud than other cluster 
galaxies (jvon der Linden et al.ll2006l ). Feedback from the dissipation of th e energy of expand 



ing r adio bubbles from such AGN has been proposed as a solution (e. g.. iBriiggen &: Kaiser 
20021 ) to the cooling-flow problem (the lack of evidence for the extreme inflows of cold gas 
which are predicted to occur in clusters due to their short radiative cooling timescales). 



Best et al.l (120061) used a sam ple of 625 nearby groups and clusters of galaxies selected 



from the SDSS (jYork et al.ll2000l ). They showed that galaxies close to the centers of clusters 
have a n enhanced probab ility of being radio-loud AGN compared to field galaxies of similar 
mass. iLin fc Mohrl (120061 ) studied a sample of 573 X-ray luminous clusters, and also found 



an enhancement of the fraction of galaxies that were radio-loud, as well as tentative evidence 
for an increase in this fraction with the mass of the cluster. 

In this paper, we use a clu ster sample more tha n 20 times larger than those used for the 
above s tudies. The rn a xBCG (IKoester et al.ll2007al ) red sequence galaxy cluster finder was 
used by lKoester et al.l (l2007bl ) to create a catalog of 13,823 galaxy clusters from the SDSS - 
the largest galaxy cluster catalog assembled to date. The clusters have redshifts between 0.1 
and 0.3, and contain between 10 and 190 red sequence galaxies brighter than 0.4 within 
the virial radius. 



We correlate the BCG positions with the FIRST (IBecker et al.lll995l ) radio source cat- 
alog, in order to study the BCGs which are radio-loud, as well as the spatial distribution of 
radio sources in the surrounding clusters. Of the 13,823 clusters in the SDSS sample, 13,240 
are within the region covered by FIRST. 



We assume an Vt m = 0.27, VL\ = 0.73, H = 71 km s 1 Mpc 1 cosmology (jSpergel et al. 



20031 ). 



2. Radio sources in the cluster environment 

We begin our analysis by finding the angular distance 8 C from the BCG position to the 
closest FIRST source. A histogram of the results for all 13,240 fields is shown in the left 
panel of Fig. [H along with a histogram of the mean number of random matches obtained 
by repeating this task around 4 positions offset by ±30' in the cardinal directions from each 
of the BCG positions. We also express the number of matches as a percentage of the total 
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number of BCGs. Note that with increasing 6 C , the BCG field counts initially fall below the 
background counts, due to the tendency of BCG fields to have a FIRST source closer to the 
BCG position than the average comparison field (which is as expected due to conservation 
of number). The background counts peak at 9 C = 147"5, which corresponds to the average 
distance from the comparison field position to the nearest FIRST source. 

The first bin (0 < 9 < 5") contains 2328 sources (compared to a mean of 7.5 ± 1.7 in 
the comparison fields); i.e., 17.5% of fields contain a FIRST source within 5" of the BCG 
position, most of which are likely associated with an AGN hosted by the BCG. An additional 
287 BCGs appear to host extended double-lobed radio sources (§EJ). Therefore at least 19.7% 
of BCGs in our sample are associated with FIRST sources. 

Some of these extended sources show up as additional counts in the 9 < 30" bins of 
Fig. [TJ In the right-hand panel, we again count FIRST sources as a function of distance, 
but this time we count all matching FIRST sources, not just the closest (so that there are 
multiple FIRST matches per field, and sometimes more than one FIRST source per annulus). 
It is easy to show that in this case the number of matches per 5" bin for the comparison 
fields is given by N = lOnprn, where p = 6.9 x 10 -6 sq. arcsec -1 is the surface density of 
FIRST sources, r is the radius of the center of the annulus, and n = 13240 is the number of 
fields. 

At 9 > 30", the counts in the BCG fields begin to fall off less steeply than at smaller 
radii, as they cease to be dominated by the intrinsic size of the central source (see §ED, an d 
begin to be dominated by other radio galaxies within the cluster. At 9 a > 400", the counts 
become indistinguishable from the background counts, suggesting that radio source activity 
is only enhanced (on average) within the central ~ 1.4 Mpc of the cluster (at the mean BCG 
redshift for this sample z = 0.2184; see also Fig. [2]). 

The power of this analysis is that even with no assumptions about the redshifts of 
sources, or attempts to associate multiple components of individual radio galaxies with the 
"parent" optical source, one can see the excess of radio sources extends far out into the 
cluster. By contrasting the left and right panels of Fig. [TJ one can also see that there is 
a significant excess of multiple-component radio sources on scales < 50", before the slope 
in the right panel changes at about 100", where the excess is mainly due to cluster radio 
sources unassociated with the BCG. 

Since the maxBCG catalog provides accurate photometric redshifts for all the BCGs, 
we can use these to convert 9 into real projected distances (and radio flux densities into 
radio luminosities; § [3]). In this analysis (Fig. [2]), we assume that each radio source in 
a particular field is at the same redshift as the BCG. As can be seen from Fig. [1] (right 
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panel), this is not true at large angular separations 9, where most of the FIRST sources 
are unassociated "background" sources. We attempt to compensate for this by scaling the 
background counts in the same manner (where we use the BCG redshift to scale the counts in 
all four associated comparison fields, and then compute the histogram of background counts 
for the entire sample). At smaller separations, most of the sources really are associated with 
the BCG, i.e., they are at the cluster redshift. It would of course be possible to match 
the other radio sources to galaxies from SDSS and attempt to select only sources at similar 
photometric redshifts to the BCG. In this analysis we wish to concentrate on the "pure" 
radio population in the clusters, however, so we defer such radio-optical matching to future 
work. 

In Fig. [21 we show the average number of FIRST sources per field, divided by the 
area of each annulus, i.e., the surface density of radio sources in the clusters. Again, there 
is a prominent peak associated with the radio-detected BCGs, which has a wide shoulder 
extending to ~ 200 kpc due to the presence of extended radio sources (see § [3]). Even out 
to large radii, there is an excess of sources in the cluster environments compared to the 
comparison fields. Since we have such a large cluster sample, we are also able to bin by other 
parameters such as radio luminosity, redshift, cluster richness, etc., and look for changes in 
the distribution of radio sources. These results will be discussed in an upcoming paper. 

In Fig . [3] we show the fraction of BCGs which are radio-loud (which, for consistency with 



Best et al 



2006 



we d efine as having an associated FIRST source with L L4 > 10 23 WHz 1 ). 
Koester et al.l (l2007bl ) provide BCG luminosities in SDSS r-band, ^-corrected to z = 0.25. 



We convert these luminosities to a bsolute magnitudes, M r , and assume an r-band stellar 
mass -to-light rat i o log jM/L) = 0.5 (IKauffmann et al.ll2003al . Fig. 14). We recover the result 
from iBest et al.l (120061 ) that the radio-loud fraction increases with increasing stellar mass. 



Best et al.l found that BCGs in rich clusters were more likely to be radio-loud than those in 
poor clusters, but concluded that this was due to the tendency of the richest clusters to host 
the most massive BCGs. Using our larger sample, we are able to show that for a given stellar 
mass the probability of a BCG to be radio loud is enhanced somewhat in richer clusters, 
suggesting that the wid e r envi ronment plays a role in radio-loudness, in agreement with the 



findings of iLin fc Mohrl fl2006l ). 



3. BCG radio luminosities and morphologies 



By converting the measured radio fluxes for the BCGs into rest-frame 1.4 GHz radio 
luminosities (-L1.4) we can investigate the dependence of radio lum inosity on the optical lumi- 
nosity of the BCG. For consistency with previous literature (e. g.. lCondonlll992l ; lLin &: Mohr 
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20061 ) we adopt an average radio spectral index a = —0.8 (S u oc z/ a ) when applying the 
fc-correction to calculate L\ A . At the relatively low redshifts considered in this paper, the 
choice of spectral index does not strongly affect our results. 



In order to account for extended emission, we use the technique of Ide Vries et al.l (120061 ) 
to find double-lobed radio sources associated with the BCGs, even where the BCG does not 
have a radio core detected by FIRST. In Fig. HJ we plot the lobe opening angl e (\1/, the angle 



subtended at the BCG position by the FIRST sources identified as lobes; see Ide Vries et al. 



for more details) for double-lobed radio sources in two size bins. In the subsequent discussion, 
we restrict our analysis to double-lobed sources with ^ > 100°, where the fraction of false 
positives (shown as the blue histogram in Fig. H]) is low. At small angles (~ 45°) we also 
see a small secondary peak in the cluster fields, due an excess (compared to the comparison 
fields) of real double sources in the cluster which are associated with galaxies other than 
the BCG. Due to their true hosts being offset from the BCG position, these are sometimes 
erroneously identified by the algorithm as being small opening-angle sources associated with 
the BCG itself (but again we emphasise that most sources with \1/ > 100° are real). We find 
453 BCGs (compared to a background expectation of 1.00) associated with double-lobed 
radio sources with \1> > 100° and diameter — 30" (83 where the core is also detected, and 
370 where it is not). Additionally, we find 196 BCGs (compared to 4.90 in the background 
case) with double-lobed radio sources with ^ > 100° and diameter 30 — 60" (108 with cores, 
and 88 without). Note that 362 of the 649 BCGs associated with doubles were also counted 
in the 2328 sources with a FIRST source within 5" discussed in § [2J 

In Fig. [5], we plot L\a (inferred from the total fluxes of the radio components associated 
with each BCG) for the 649 double-lobed sources (red squares) as a function of M r . The 
black squares in Fig. [5] are the 1966 BCGs not identified as doubles, which nevertheless have 
a FIRST sour ce within 5". We a lso pl ot the empirical division between FR I and FR II radio 
galaxies from iLedlow fc Owenl (1l996[ ) , where we have assumed M r = M24.5 (the isophotal 
magnitude to 24.5 mag arcsec -2 from lLedlow fc Owenl ). 



In contrast to the results of Ide Vries et al.l (120061 ) for SDSS quasars, we find very little 
excess at ~ 180° for diamet ers 60 — 120" ( i.e., t here are very few doubles as large as 
this in our BCG sample). Since Ide Vries et al.l (120061 ) studied quasars with median redshift 
z = 1.40 and associated radio emission from FIRST, their sample has a higher mean radio 
luminosity, so this result is perhaps unsurprising, even given the smaller angular diameter 
distance to ou r sources; the p hysical sizes of our sources are about a factor 4 smaller than 
the FR lis of Ide Vries et all peaking at 55kpc (Fig. E]), and dropping to zero at around 
200 kpc. As can be seen from Fig. (2J 200 kpc is where the slope of the histogram changes, 
showing that the excess matches at > 200 kpc are dominated by other radio galaxies in the 
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cluster, and at < 200 kpc, by the presence of multiple-component radio sources associated 
with the BCG. 

Fig.[5]shows that most o f the double-lobed sources are located near the FR I / II division 
from iLedlow fc Owen! (119961 ) , which adds weight to the argument that radio AGN in cluster 
environments at low redshift do not tend to be powerful FR lis ( jHill fc Lilly! Il99ll ) - they 
could certainly be seen by us if present. It is clear that the upper envelope in radio luminosity 
in this plot shows an increase in radio luminosit y with optical lumi nosity (a proxy for stellar 



mass), in agreement with previous results (e.g., iJarvis et al.ll200ll ). 



We note that only a few detected sources (7.2%) are fainter than 10 23 W Hz -1 (essentially 
due to the FIRST flux limit). For s uch a radio luminos ity to be due entirely to star formation 
would require SFR ~ 100M o yr _1 (IHodge et al.ll2007l ). so we conclude that radio- loud AGN 
are the dominant source of radio emission in those of our BCGs detected in FIRST. Of course, 
even in sources with S\a ^> 1 mJy, star formation may contribute to the radio emission. In 
some cases, this star formation will also manifest as optical emission lines, whereas in other 
cases the star formation may be optically obscured. 



4. Summary and future work 

By matching the maxBCG and FIRST catalogs, we were able to determine that 19.7% 
of BCGs in our sample are detected in the radio; the majority (92.8%) of these are radio- 
loud AGN (L1.4 > 10 23 WHz -1 ). The radio loud fraction has a strong dependence on stellar 
mass of the host, increasing from ~ 1% at 4 x 1O 1O M to ~ 30% at 5 x 1O 11 M . For a 
given stellar mass, BCGs in rich clusters have a higher tendency to be radio-loud than BCGs 
in poor clusters or groups, suggesting that the wider environment has some influence on 
radio-loudness. 

Signs of double-lobed emission are seen (at the resolution of FIRST) in 24.8% of the 
2615 BCGs associated with FIRST sources. The physical size of these double-lobed sources 
peaks at 55 kpc, with a tail extending to ~ 200 kpc. Their radio luminosities place them 
close to the FR I / II division. 

There is an excess of radio sources in our clusters, which are associated with other cluster 
members, extending out to ~ 1.4 Mpc, but the central BCGs are, overall, the dominant source 
of AGN activity in these clusters. 

Future investigations using this dataset could address the relationship between optical 
spectroscopic indicators of star formation and AGN activity (as determined from the SDSS 
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"value added" catalogs of iBrinchmann et al.l 12004 and iKauffmann et al.ll2003bl ). radio loud- 
ness / luminosity, and environment, as well as extending the FIRST-SDSS match to other 
galaxies in the clusters (as determined by their photometric redshifts). 



For those BCGs unde tected in FIRST, it is possible to "stack" the FIRST maps (IWhite et al 



20071 ; Ide Vries et al.l 120071 ) to determine the average radio properties of sources too faint to 
be detected individually. In this way, we can study BCGs with radio emission two orders of 
magnitude below the FIRST flux limit. These results will be discussed in a future paper. 



Work was performed under the auspices of the U. S. Department of Energy, National 
Nuclear Security Administration by the University of California, Lawrence Livermore Na- 
tional Laboratory under contract No. W-7405-Eng-48. Research made use of the VLA FIRST 
Survey (NRAO programs AB628, AB879 and AB950) and the SDSS. The National Radio 
Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 
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Fig. 1. — Histogram of the total number of FIRST sources in all 13,240 BCG fields, as a 
function of angular distance 9 (in arcseconds) from the target (BCG or offset comparison 
field) position, in bins (i. e., annuli) of width 5". The black histogram shows the counts in the 
fields containing the BCG, and the grey histogram shows the mean counts in the comparison 
fields. The left panel shows just the closest FIRST source to the target position (i.e., for 
each field, we consider only the single FIRST source closest to the target), whereas the right 
panel shows all FIRST sources in the target field (which can occasionally include multiple 
FIRST sources in a given annulus). The right-hand vertical axis in both panels shows the 
counts per bin as a percentage of the total number of fields (13,240). 
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Fig. 2. — Histogram of the average projected surface density of FIRST sources (black), as a 
function of distance in kpc from the BCG position (assuming, for each field, that all FIRST 
sources are at the BCG photometric redshift). The grey histogram shows the comparison 
field counts, at the mean of four positions offset by 30' from the BCG positions. 
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Fig. 3. — The percentage of BCGs that are radio-loud AGN, as a function of stellar mass 
(black line). Also plotted are the same data split into two subsamples; BCGs in poor clusters 
(< 20 members; blue line) and rich clusters (> 20 members; red line). At a given stellar 
mass, BCGs in rich clusters have a slightly higher probability of being radio-loud. 




Fig. 4. — Histogram of opening angles for double and triple radio sources associated with 
BCGs (black) and expected contamination from random superpositions of sources which 
appear to be doubles / triples (grey). The left panel shows doubles / triples with major axes 
between and 30", and the right panel, major axes between 30 and 60". The bin size is 5°. 
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Fig. 5. — Rest-frame 1.4 GHz radio luminosity as a function of SDSS r-band absolute mag- 
nitude for the BCGs detected in FIRST, fc-corrected to z = 0.25. BCGs with a FIRST 
source within 5" are shown as black squares, except for cases where our "double-detector" 
algorithm (see § [3]) finds a double radio source associated with the BCG (sometimes with 
a radio-detected core, sometimes without). These 649 doubles and triples are plotted as 
red squares, with the radio luminosity determined from the total integrated radio flux of 
the components. The s olid line is the division between FR I and FR II radio galaxies from 



Ledlow fc Owen! (119961 ) (we make no correction from their M24.5 to our M r 
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Fig. 6. — N d , the number of double- lobed radio sources as a function of physical size of the 
radio emission for the 649 doubles in our sample. 



